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ABSTRACT 

We examine the effect of inhomogeneous reionization on the galaxy power spectrum and the 
consequences for probing dark energy. To model feedback during reionization, we apply an 
ansatz setting the galaxy overdensity proportional to the underlying ionization field. Thus, in- 
homogeneous reionization may leave an imprint in the galaxy power spectrum. We evolve this 
imprint to low redshift and use the Fisher-matrix formalism to assess the effect on parameter 
estimation. We show that a combination of low- (z — 0.3) and high- (z = 3) redshift galaxy 
surveys can constrain the size of cosmological HII regions during reionization. This imprint 
can also cause confusion when using baryon oscillations or other features of the galaxy power 
spectrum to probe the dark energy. We show that when bubbles are large, and hence detectable, 
our ability to constrain w can be degraded by up to 50%. When bubbles are small, the imprint 
has little or no effect on measuring dark-energy parameters. 

Key words: cosmology: theory - cosmological parameters - galaxies: formation - inter- 
galactic medium 



1 INTRODUCTION 

During the epoch of reionization, groups of star-forming 
regions generate significant numbers of ionizing pho- 
tons, which may lead to HII regi ons many Mpc in size 
JFurlanetto. Zaldarriaga & Hernq uist 2004). These ionized bubbles 
grow as further structure forms, eventually merging and causing 
full reionization of the intergalactic medium (IGM). Conditions 
within these HII regions may be significantly different than in the 
surrounding neutral IGM. For example, the temperature in these 
HII regions will be raised by photoionization heatin g, which is 
known to suppre s s star formation in low-ma s s haloes l Ree sill 9861: 
Efstathioul 1 19921: iThoul & Weinberg! 1 19961 : iKitavama & Ikeuchil 
2000; Dii kstra et alj |2004). Also, the ionizing flux generates 
more free electrons, which af fects the abundance of molecular 
hydrogen JOh & Haimanll2002l) . an important coolant. These, and 
other feedback mechanisms, will affect the fraction of baryons 
that condense in haloes, and in turn modify the number density 
of directly observable galaxies dBarkana & Lo eb 2001). This 
suppression will be inhere ntly inhomogeneous, as highly biased 
regions will ionize first JBabich & Loebl 120061) . Understanding 
the detailed effects of feedback is one of the major remaining 
challenges in understanding galaxy formation. 

Galaxies formed during reionization will be low in mass and 
faint by comparison to those from later generations of galaxy for- 
mation. This, along with absorption along the line of sight by 
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the IGM, makes it difficult for existing telescope facilities to de- 
tect large numbers of early galaxies directly. However, a num- 
ber of large galaxy surveys now exist, which probe the distri- 
butio n of galaxies in th e lower-redshift regime I IEfstathiou et alj 
120021 : ISeliak et alj 120051) . These surveys focus on high-mass lu- 
minou s objects - e.g ., luminous red galaxies (LRGs) in SDSS 
(Eisenstein et al. 2005) - which are easily detected. The abundance 
of such objects will depend in a non-trivial way upon the number 
of low-mass progenitors, especially upon the amount of condensed 
gas available for mergers. Thus, these late-forming galaxies will, 
indirectly, be affected by the efficiency of galaxy formation during 
reionization. Motivated by these arguments, we consider the possi- 
bility that large galaxy surveys in the low-redshift Universe may be 
used to probe inhomogeneous reionization through its feedback on 
early galaxy formation. The many uncertainties remaining in our 
understanding of galaxy formation make it difficult to develop a 
rigourous formalism for this imprint, and motivate a simpler, hope- 
fully more robust, approach. 

Besides the possibility of detecting reionization, its imprint 
in the galaxy power spectrum may act as a source of noise when 
probing cosmology. M odern observations of the c o smic microwave 
background (CMB ) (de Bernardis et al. 2000; Halvers onet all 
l2002l;lMason et al . 2003; Benoit et al. 2003; Goldstei net"ai1l2003l : 



etaLp 
20031) " 



Spergel et al. 2003) have greatly extended our knowledge of cos- 
mological parameters. One result has been the realisation that 
~ 70% of the Universe is composed of an unknown form of en- 
ergy that g enerates the acce l erated expansion seen in SN la ob- 
servations JRiess et alj|l99a : IPerlmutter et alj|l999l) . This is one 
of the most puzzling discoveries of our times, and it is hoped 
that future observations in the fields of SN la iRiess et alj|l998l 
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l2004|Perlmutter et al.ll 1999b. weak lensing JHoekstra et alj|2005b. 



and galaxy surveys l lSeo & Eisen stein 20031; iBlake & Glazebrookl 
2003) will constrain the time evolution of the dark energy giving 
clues as to its nature. For this reason, in this paper, we will focus on 
how reionization may affect estimates of dark-energy parameters. 

Large galaxy surveys contribute information on the dark 
energy in two main ways. First the form of the matter power 
spectrum, probed by galaxy surveys via the proxy of the 
galaxy power spectrum, depends on different parameter com- 
binations t han the CMB, breaking many o f the parameter de- 
generacies jEisenstein, Hu & Tegmarklll999b . Secondly, the pre- 
recombination oscillation of the photon-baryon fluid leaves an 
imprint in the matter power spectrum, which may be used as 
a standard ruler to determine the an gular diameter d i stance 
Da(z) as a function of redshift z JSeo & Eisensteinl 1200 3l ; 
IBlake & Glazebrook 2003). T hese baryon oscillatio ns have now 
been detected fcole et alteOOStlETsenstein et alj2005b by both 2dF 
and SDSS. If the imprint in the galaxy power spectrum from patchy 
reionization can mimic or conceal any feature of the galaxy power 
spectrum from density fluctuations, then our ability to constrain 
dark energy using galaxy surveys will be degraded. 

In this paper, we explore the possible consequences of this 
environmental dependence on the galaxy power spectrum. The 
process of galaxy formation is still only poorly understood and 
so a detailed analysis of feedback is inappropriate. Instead we 
choose to link galaxy formation to the neutral fraction by a sim- 
ple ansatz, by which we hope to bring out the underlying be- 
haviour, leaving the details for a later age. In keeping with this 
"simple is best" ideology, we choose to model the variation in neu- 
tral fr action using an analogue of the halo model dCoorav & Shethl 
2002). With this approach we hope to phrase the problem in a 
general fashion, avoiding detailed assumptions about the reioniza- 
tion history. To address these questions i n a quantitative fashio n, 
we employ the Fisher-matrix form alism dJungman et al.l l 1996a b; 
iTegmark. Taylor & Heavensll 1997b . This allows us to convert a the- 
oretical model into predictions for the parameter constraints attain- 
able by imagined experiments. 

The outline of this paper is as follows. In Sj2]we detail the form 
of the ionization power spectrum and describe our simple ansatz re- 
lating it to galaxy formation. Then, in Sj3] we bring the two together 
describing the complete model galaxy power spectrum, including 
the effects of redshift distortions and the Alcock-Paczynski effect. 
In !j4] we outline the Fisher-matrix formalism. Having set out our 
model, in ij5]we discuss the possibility of detecting reionization us- 
ing galaxy surveys. This is then expanded to consider the implica- 
tions for dark-energy constraints in SJ6] Finally, in iJ7]we summarise 
our conclusions. 



2 BUBBLE MODEL 

We wish to relate the overdensity of galaxies to the ionization frac- 
tion within a given region. To do this, we make the simple ansatz 
that there is a component to the galaxy power spectrum which lin- 
early traces the ionized fraction. Hence, we may write the number 
density of galaxies n(r) at position r as 



n(r) = n[l + 5 ga a(r) + <5 bub (r)], 



(1) 



where n is the mean number density of galaxies, 5 ga i(r) = b<5(r) 
assumes galaxies trace the underlying dark-matter fluctuations 8 
with bias b, and we calculate the fractional overdensity of galaxies 



due to an ionization field x%(r) by (5bub(r) = —ebXi(r), where et 
parametrizes the strength of the effect. 

Writing the number density in the form of equation ([TJ leads 
to a galaxy power spectrum 



P(k) 



(l-e b Q)' 



[Pgal(fc) + 2P ga i, bub (fc) + P bub (fc)] , (2) 



where Q is the filling fraction of the bubbles. For simplicity, we 
choose to neglect the cross-correlation, which will be smaller or 
comparable in size to the other terms and represents an unnecessary 
refinement given the simplicity of our toy model. Note the overall 
rescaling of the power spectrum because the mean galaxy density 
(n) = n(l — ebQ). Typically, e b Q <C 1 and we can neglect this 
correction and take 



P(fc) = P ga ,(fc) + P bub (fc). 



(3) 



We now need to calculate the bubble power spectrum 
Pbub(fc). In order to phrase the problem as broadly as possible, 
we eschew detailed assumptions about reionization in favour of a 
more general approach. In this paper, we choose to associate re- 
gions of ionization with "bubbles," in analogous fashion to the halo 
model's as sociation of mass with h aloes. Following the halo-model 
formalism l lCoorav & Sheth 2002), P bub (fc) is given by the sum of 
two terms, P bub (fc) = P lb (k) + P 2b (k), which describe corre- 
lations within the same bubble and between two different bubbles 
respectively. These terms are given by 

2 



P lb {k) 



t b J dmn(m) ( — ) \u(k\m)\ 2 , (4) 



P (k) = t h I dm\ n(mi) I —=- ) u{k\m\) 

x / dm2n(m2) f —=- 1 u{k\m-i)Pbb{k\'mt,mi), (5) 

where n(m) is the comoving number density of bubbles of mass 
m, Pbb(k\mi, 7712) is the power spectrum of bubbles of mass mi 
and mi, and u(k\m) is the Fourier transform of the bubble ion- 
ization profile u(r\m). With this notation, we may write the vol- 
ume filling factor of the bubbles as Q — J dmn(m)(m/p), and 
the bubble volume as Vbut> = m/p — 47rr bub /3, where r bub is 
the comoving bubble radius. Throughout this paper, we will as- 
sume a top-hat profile u(r\m) = 0(|r| — r bub ))/V bub , for which 
u(k\m) = 3ji (krhub) / (krhub) , where jt(x) is a spherical Bessel 
function of order I. 

If we assume a delta-function size distribution and treat the 
power spectrum of the bubbles as tracing the dark-matter power 
spectrum Pbb(k) « P&s(k,z — z r i), where z r i is the redshift at 
which the imprint is formed, this reduces to 



P 16 (fe) 



e b QVbub\u(k\m)\ 



P 2b (k) = e 2 b Q 2 \u(k\m)\ 2 P ss (k\m). 



(6) 



(7) 



In order to keep our model simple, we ignore evolution in the 
bubble-size distribution. In reality, the relevant bubble sizes will 
be determined by the period when most baryons condense, an 
extended process that will average over the evolution of bubble 
growth. We also ignore the effects of bubble overlap, which is ex- 
pected to occur for large Q and undermines the halo-model ap- 
proach. Once bubbles begin to overlap, using isolated spheres to 
model the HII regions will not correctly represent the true size 
and shape of the ionized regions. To a first approximation though, 
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Figure 1. Comparison of the galaxy and bubble power spectra. We plot 
Pgal(ft) at two redshifts z = 0.3 (solid curve) and z = 3.0 (long 
dashed curve). For each redshift, we plot the non-linear scale: fc max (2 = 
0.3) = 0.11ft Mpc" 1 and fc max (z = 3.0) = 0.53/iMpc- 1 (dashed 
vertical lines, from left to right). For comparison, we plot Pbub(^) f° r 
the parameters (rbub = 80Mpc,£b = 0.6) (short dashed curve) and 
( r bub = 20Mpc,tb = 0-6) (dotted curve). Notice how the latter curve 
resembles constant white noise in the region k < O.l/iMpc^ 1 . The for- 
mer curve displays a cutoff in power close to the galaxy-power-spectrum 
peak. Finally, we plot a bubble power spectrum Pfzh(^) (dot-dashed 
curve) that has been calculated using a bubble size distribution taken from 
iFurlanetto eta?] J2004I) . with (r bub ) ps 20 Mpc. 



this effect will give an effective distribution of bubble sizes, and so 
should not affect our qualitative conclusions. We will take Q — 0.5 
and z r i — 6 in what follows, and use (rbub, £(>) to parametrize the 
bubble power spectrum. Note that the 2-bubble term is subdomi- 
nant in the regime that we consider, making the details of z r i, and 
any biasing of Pbb(k) with respect to Pgs(k), unimportant. 

Figure Q] shows the form of the bubble power spectrum in 
this model. Note that the power is fairly constant on small k and 
cuts off sharply on linear scales smaller than the bubble radius. 
As a simple example of including a smooth bubble size distri- 
bution, Figure Q] shows P bub c alculated using the bubble distri- 
bution of Furlan etto et alj q2004), assuming an ionizing efficiency 
C, = 40 and Q = 0.83, which gives a volume averaged bubble size 
(rbub) ~ 20 Mpc. The main effects of the distribution in bubble 
sizes are to smooth out the oscillations seen in the single size model 
and to decrease the rate at which power decreases on scales below 
the characteristic bubble size. Having shown that the high-fc cut- 
off occurs even with a smooth distribution of bubble sizes, we will 
henceforth restrict ourselves to the simpler, single bubble size case. 

The onset of non-linearity limits the scales that a galaxy sur- 
vey is able to probe. We choose to define this cutoff scale by requir- 
ing that the average fluctuation on a scale R satisfies cr(R) < 0.5 
for R — 7r/(2fc max ) dSeo & Eisensteinll2003r) . This cutoff can lead 
to a degeneracy between the bubble spectrum and the constant shot- 
noise expected on large scales from non-Gaussian clustering of the 
galaxies (Seljak 2000). We see in Figure[Tjthat on large scales the 



bubble power spectrum becomes constant. If, for a given galaxy 
survey, rbub is sufficiently small, then the curvature of the bubble 
spectrum will lie at k > fc max , and the bubble spectrum will be in- 
distinguishable from shot-noise. Including galaxy surveys at higher 
Z, where the non-linear scale is smaller, helps break this degener- 
acy. 

We note that, for a random variable with zero mean, we would 
expect the power spectrum to vanish on large scales. That this does 
not occur relates to a generic problem of the halo model 1-halo 
term, which is constant on large scales. In most applications this 
is masked by a dominant 2-halo term, which does decrease on 
large scales. However, in our model the 2-bubble term is negligible 
making this issue obvious. Given that our model predicts a bub- 
ble power spectrum that looks like shot-noise on large scales, we 
must worry both about how the removal of shot-noise will affect 
our results and how to distinguish the effect of bubbles from shot- 
noise. As mentioned above, the existence of a cutoff in the bub- 
ble power spectrum distinguishes it from shot-noise (although we 
must observe this cutoff for this to work). When we come to anal- 
yse the effect of bubbles on cosmological parameter estimation, we 
will include a term representing white shot-noise to account for this 
possible confusion. 

Having generated an imprint from patchy reionization, we 
must evolve it to lower redshift. Our knowledge of how mergers 
recycle matter from many smaller haloes into fewer more massive 
haloes is not sufficient to handle this rigorously. Instead, we will 
consider three cases that ought to bracket the truth. We take rbub to 
be a constant, fixing the shape of Pbub(k), and then consider how 
its amplitude varies with time. We will consider three models for 
this time evolution 



PbubV"'? z — Zri), 
(U.ik.:) = { *WM = 0) [^] 

P b ub(M = ^) [- eki 



G(z = z ri ) 



(Model A), 
(Model B), 

(Model C). 



In Model A, we assume that, once produced, the power spectrum 
Pb\ib{k) from bubbles remains constant in time. As the density 
fluctuations continue to grow this means that Pbub becomes less 
significant at later times. In Model B, we allow Pbub(fc) to grow 
as the square of the linear growth function G(z). Thus in Model 
B, Pbub(k) remains a constant fraction of the total galaxy power 
spectrum. It seems unlikely that the bubble imprint would grow in 
this fashion, but we include this model in order to consider the case 
where the bubble imprint is equally important at all redshifts. Note 
that in this model, we choose to normalise the bubble spectrum to 
the present day. This provides a simple way of restricting Pbub(fc) 
to amplitudes comparable to the density power spectrum. Finally, 
with Model C, we consider the case where Pbub(fc) decreases with 
time. This will provide an estimate of the worst-case scenario for 
detecting the bubbles. The time evolution of Pbub(fc) is most im- 
portant when we can compare surveys at different redshifts. In the 
case of a single redshift survey, any growth can be absorbed into an 
effective et, for that survey. 

What range of values may our two free parameters et 
and rbub reasonably take? The charact eristic size of the bub- 
bles is the easiest question to address. IFurlanetto et alj (2004) 
present a model for bubbles forming around highly biased re- 
gions lea ding to typical sizes of ~ 5 Mpc, wh en Q = 0.5 
(see also Furlanetto, McOuinn & Hernquistl I2006T) . In contrast, 
IWvithe & L oeb (2004) use arguments based on light-travel times 
and cosmic variance to obtain bubble sizes of ~ 60 Mpc at the end 
of reionization. This latter value can be taken as an upper limit on 
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reasonable bubble sizes, while the former gives a more reasonable 
estimate of what we might expect. These values a re in broad agree- 
ment with the r esults of computer simulation JHiev et al.l 120051 : 
IZahn et alj2006l) . which yield sizes ~ 10 Mpc. 

The range of eb begs the question of how exactly to inter- 
pret this parameter. We have assumed a linear relation between the 
ionization fraction of a region and its galaxy overdensity. We can 
readily see that n(r) > 0, which implies a solid upper limit of 
£(, < 1. An alternative approach is to consider the suppression of 
galaxy forma tion in haloes of low mass. Simulations at low red- 
shift (z < 3) JThoul & Weinberg! 19961 : 1 Kitavama & IkeuchfeOOOh 
indicate significant suppression of galaxy formation in haloes with 
circular velocities Kirc < 50kms~ . At higher z, photoioniza- 
tion is less effective due to the decreased cooling time, decreased 
UV flux, increased self-shielding from the higher densities, and 
collapse beginning b efore any UV background can be generated 
JDiikstra et ai]|2004l) . In this case, iDiikstra et all J2004I) find that 
only haloes with V c i rc < 20kms _1 suffer significantly reduced 
condensation of baryons. To estimate the mass fraction in galaxies 
affected by photoionization feedback, we take this latter value and 
apply it to the Press-Schecter distribution ( Press & Schechter 1974) 
as a low-mass cutoff below which no galaxies form. This gives an 
estimate of the decrement in galaxies due to photoionization feed- 
back, 



66 



A 



F(M > Affcedback) ~ F(M > Mcool) 

F(M > Mcool) 



(8) 



where F(M) is the fraction of mass in haloes of mass greater than 
M, Mfeedback is the mass corresponding to Vcirc = 20kms -1 , 
and M coo i is the mass corresponding to the virial temperature 
T v i r ~ 10 4 K needed for effective cooling by atomic hydrogen. 
Evaluating equation l[8} gives ej, ft* 0.18 at z — 10 and et ~ 0.10 at 
z — 6, which give an indication of sensible values. Once a galaxy 
grows large enough, gravity will overcome feedback of this form 
and damp this effect. Thus, these numbers represent an effective 
upper limit in the most plausible model. 



3 GALAXY POWER SPECTRUM 



In 



constructing our g alaxy power spectrum, we follow 



& 

ISeo & Eisensteinl J2003D . Incorporatin g the effects of bias, 
linear redshift distortions ( IKaiseg|l987h . and linear growth, the 
galaxy power spectrum takes the form 



Pgal(k, fj,,z) = 



G(z) 



G(z = 0) 



b 2 (l + Pn 2 ) 2 P s {k,z = 0), (9) 



where P$(k, z = 0) is the power spectrum of the dark matter at 
the present day, b = f2 m (z)°' 6 //3 is the bias, and /i 2 = k 2 Jk 2 
is the direction cosine between the Fourier-mode wavenumber and 
the line of sight. We define the redshift-distortion parameter f3 in 
terms of as, 9 and as, the fluctuations in galaxies and dark matter, 
respectively, smoothed on scales of 8k" 1 Mpc, by the relation, 



<x 8 , 9 = agb^/1 + 2/3/3 + /3 2 /5. 



(10) 



In order to calculate the linear growth factor G(z), we integrate the 
perturbation equation, 



G(t) + 2HG(t) - 4nGp m G(t) = 0, 



(11) 



with 
H 



-2 = o»(i + *r + (i - n™ - n x )(i + z y + n x ( z ), (12) 

tin 



and where the energy density in dark energy is given by 



fix (z) = fix exp 



3 / dz 



, l + w(z') 



1 



(13) 



In the special case of a cosmological constant, the growth factor 
may be expressed as 



g(z) = -a 



H{z) 
Ho 



1 + z' 



[H(z')/H ]3 



dz , 



(14) 



but for a general dark-energy model where w(z) ^ — 1, the 
full numerical integration is nece ssary JWang & Steinhardtll 19981 : 
IWeinberg & Kamionkowsklll2003l) . 

Figure Q] shows P ga i(fe,/i, z) averaged over angle and eval- 
uated at z — 0.3 and z = 3. It displays a clear peak at 
k ftJ 0.02/tMpc -1 , corresponding to the scale of matter-radiation 
equality, and visible baryon oscillations on smaller scales. These 
features arise from the acoustic oscillation of the baryon-photon 
fluid during the period of tight coupling before recombination. The 
sound speed, which governs the peak positions, is well measured 
from the CMB. Consequently, the baryon oscillations may be used 
as a standard ruler, allowing a direct measurement of the angular- 
diameter dhtanc^^he^efe^turesjiavenow been detected in galaxy 
surveys dCole et al.ll2005l : Eisenstein et al. 120051). and their use in 
probing the dark energy i s well known l lSeo & Eisensteinl 120031 
iBlake & Glazebrookll2003l) . 

When converting the observed redshift and angular position 
of galaxies into linear space, we must assume a particular cos- 
mology. If this reference cosmology is different from the true cos- 
mology, then we will introduce distortions into the inferred dis- 
tribution of galaxies. This is the Alcock-Paczynski (AP) effect 
I IAlcock & Paczvn ski 1979) and is essentially a cosmological red- 
shift distortion. We may express the power spectrum inferred by 
our observations in terms of the true power spectrum P tI (k tI , pi tr ) 
by 



Pobs(k, (j) 



>a\z 



)H % 



D t X 2 {z)H{z 



-P^fc'V 1 ), 



(15) 



where H and Da are calculated using the reference cosmology, and 
H tT and D X X with the true cosmology. We write the components of 
a Fourier wavevector parallel and perpendicular to the line of sight 
as k\\ = {H^/H)^ and k± = (L>a/£>a>-L- The information 
contained in the AP effect can be useful in probing the evolution of 
the dark energy, so we include it in this analysis. 

The final observed galaxy power spectrum incorporates all of 
the effects that we have discussed before and takes the form 



P obs (k, n) = -f 



D\{z)H^{z) 



DT^)H(z) 



P gal (fc tr , M tr )+P bub (fe t 



+ Pshot, (16) 



where Ps hot is residua l shot-noise from non-Gaussian clustering of 
galaxies (Seljak 2000), which we treat as a constant white-noise 
term. 



4 FISHER MATRIX 

To quantitatively constrain the effect of bubbles on the galaxy 
power spectrum, we turn to the Fisher matrix. This formalism al- 
lows us to estimate the uncertainties on a set of model parameters 
= {61,62, ...,6n) given some data set. We define the Fisher 
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matrix JTegmark et alJI 19971) 



± 1 1 



2 log L 
09,09, 



(17) 



where L is the likelihood function describing the probability dis- 
tribution of the parameters and Oo is the place in parameter space 
where the Fisher matrix is evaluated, typically the point of maxi- 
mum likelihood. Given the Fisher matrix, the Cramer-Rao inequal- 
ity states that the minimum uncertainty on a parameter 9i is given 
by A9i > (F~ \[ . This estimate of the uncertainty will be reli- 
able provided that Oo is near to the true values of the parameters. 

To evaluate Fij , we need to specify a model, which determines 
the dependence of the likelihood function on O, and a point in pa- 
rameter space where we wish to determine parameter uncertainties. 
In the case that the model parameters are Gaussian distributed, the 
Fisher matrix takes the form 



F 



itr(C- 1 C, a C- 1 C, /3 ) + ^^( l 



0/j, 



09^ 



(18) 



where C is the covariance matrix for the data, and /i is the data's 
mean. This will be a good approximation in the case of both CMB 
observations and galaxy survey. Note that, for our purposes, we 
will need to combine information from both the CMB and galaxy 
surveys. When used together these data sets break many degenera- 
cies that are present when they are used alone. Let us consider the 
Fisher matrix from each of these in turn. 

A CMB experiment may be characterised by a beam 
size #bcam and sensitivities to temperature ot and polariza- 
tio n a p. Given these quantiti es, the Fisher matrix is given 
by ( Jungman et al. 1996a b; Kam ionkowski. Kosowskv & Stebbinj 



Zaldarriaga & Seliakll 19971) 



EE^gWi 



ac) 



(19) 



where Cf is the power in the fth multipole for X = T, E, B, 
and C, the temperature, E-mode polarization, B-mode polariza- 
tion, and TE cross-correlation respectively. The elements of the 
covariance matrix Covi between the various power spectra are 
JKamionkowski et al.lll997HZaldarriaga & Seliaklll997h 



(Cov«)tt = 

(Cov^)ee = 

(CoVf) B B = 

(Cov<?)cC = 

(Covi)te = 

(Cov^)tc = 

(Covf) E c = 

(Cov^)tb = 



(Cti +w T 1 B e 2 ) 2 , 



(2£ + l)/ sky 
2 



(C E i + w P B e ) , 



(2£ + l)/ sky 

(2£ + l)/ sky ( ^ + ^ 1 ^ 2)2 

[C C i + (Cti + w T B\ 



(2£ + l)f sky 

re 
2 



xiCTt+w^Bz 2 )], 



Ccii 



CciipTi + w T B e ), 



(2£ + l)/ sky 

2 
(2£ + l)/ sky ' 

(op.^f — C C i(C E e + w P 1 Be 2 ), 

(Cov^Jeb = (Cov,.)cb = 0. 



(20) 



Here B t is the beam window function, assumed Gaussian with 
Bj = exp[-£(£ + l)0b eam /8 In 2], where 6> bcam is the full-width, 
half-maximum (FWHM) of the beam in radians. Also, wt and wp 



Table 1. Specification for CMB experiments 



Experiment 


Frequency 


"beam 


(Tt; 


Op 


WMAP 


40 


28.2 


17.2 


24.4 




60 


21.0 


30.0 


42.6 




90 


12.6 


49.9 


70.7 


Planck 


143 


8.0 


5.2 


10.8 




217 


5.5 


11.7 


24.3 



Notes: Frequencies are in GHz. Beam size 6>t,cam is the FWHM in arcsec. 
Sensitivitie s ctt and op are in /xK per FWHM beam, w = (0t, cam cr) — 2 ■ 
Taken from Eisens tein et alj U999I) . 

are the inverse square of the detector noise for temperature and po- 
larization, respectively. For multiple frequency channels we replace 
wtB 2 with the sum of this quantity for each of the channels. 



F, 



Moving no w to galaxy surv eys, we may write the appropriate 
markfl997l) 

<9mP(k)dlnP(k) 



axy 
Fisher matrix as CTegmarkl l 19971) 



-K ff (k) 



d J k 



(21) 



lo 09 a 09,3 y '2(271)3 

where the derivatives are evaluated using the cosmological param- 
eters of the fiducial model and V a s is the effective volume of the 
survey, given by 



V e s{k,n) 



n{r)P{k,n) 

n(r)P(k,fi) + 1 



nP(k,/j,) 



nP(k,n) + i 



v„ 



(22) 



Here the galaxy survey is parametrized by the survey volume 
Kurve y and the galaxy density n(r), which in the last equality 
we assume to be uniform n. In addition, we must specify /c max , 
a cutoff on small scales to avoid the effects of non-linearity. We 
choose to define this cutoff scale by the criterion a(R) < 0.5 for 
R = 7r/(2fc max ) JSeo & Eisensteirj|2003h . 

To apply the above framework, we need a theory relating the 
observables C; X and P(k, /i) to the parameters. For the CMB this 
is standard, while in the case of the galaxy surveys we use equation 
J16t , which arose from our discussion in Sj2]and i|3] Using these 
models, we calculate the Fisher matrices for individual galaxy sur- 
veys and our CMB experiment, and then combine them 



7 a0 



F 



a/3 



+E^. 



a/3 



(23) 



where i labels the different galaxy surveys. This total Fisher matrix 
is then inverted to get parameter error predictions. 

For this calculation, we need the specifications of our experi- 
ments. These are given in Table s[TJ and [2] We consider two galaxy 
surveys. The first uses parameters corresponding to the SDSS LRG 
survey, which is currently underway. The second is a hypothetical 
survey at z = 3, based upon a survey of Lyman break galaxies 
dSeo & Eisensteinll2003l) . 

Finally, we must decide upon our choice of parameter space. 
We specify our cosmology using 7 parameters describing the total 
matter fraction Q, m h 2 , Q, m , the baryon fraction Q,bh 2 , the inflation- 
ary amplitude A 2 S , scalar spectral index n 3 , optical depth to last 
scattering r, and the tensor-scalar ratio T/S; each galaxy survey is 
described by five parameters (log H, log Da, log G, log j3, P a hot)\ 
to these, we add two parameters (et, rb u b) to describe our bub- 
ble model (we choose to set Q = 0.5 and z r j = 6). I n cho osing 
these parameters, we are following ISeo & Eisenstein! J2003I) . We 
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Table 2. Specification for galaxy surveys 



Survey 


z 


•'survey 

0~ 3 Gpc) 3 


n 

(h 3 M P c~ 3 ) 


k 
™max 

(/iMpc" 1 ) 


°"8, 9 


SDSS 


0.3 


1.0 


10" 4 


0.11 


1.8 


S2 


3.0 


0.50 


10~ 3 


0.53 


1.0 



Notes: Taken from Seo & Eisen stein 12003). 
1 




20 40 60 80 100 

W m p c ) 

Figure 2. Model A: Contour map of detection e^ > 2cr eb and r^ub > 
2oy bub in the bubble parameter plane. The white region is detectable by 
SDSS alone, the grey region is detectable by SDSS+S2, and the black re- 
gion is undetectable to all surveys. Planck CMB data is assumed in all cal- 
culations. 



treat all of the above parameters as being independent and then ex- 
tract information about the dark energy from our uncertainties on 
(log H, log Da) from each survey. We choose to parametrize the 
dark energy using three parameters (Six, Wo.Wi), taking the dark- 
energy equation-of-state parameter to be w(z) = wo + w\z. In 
deciding on our fidu cial values, we follow the results of WMAP 
dSpergel et alJl2003T) for the cosmological parameters. These are 
broadly consistent with the updated results of Spergel et al. ( 2006), 
except for the decreased J2 m and t. Evaluating the Fisher matrix 
at these different best-fitting parameters modifies our constraints 
only slightly. The bubble parameters are highly uncertain, and so 
we choose to explore a large parameter space. 



5 POSSIBILITY OF DETECTING BUBBLES 

Now that we have established a theoretical framework, we wish to 
determine whether the imprint can be detected using the specified 
galaxy surveys. Our null hypothesis is that there is no imprint, and 
we assume that a detection requires that we can distinguish both e;, 
and rbub from zero at approximately the 2<j level; i.e., we require 
both eb > 2<T ei) and rbub > 2oy bub . Throughout, we assume the 
inclusion of CMB information at the level of Planck. Less precise 
CMB data will relax constraints on cosmological parameters, caus- 
ing parameter degeneracies to decrease the sensitivity of the galaxy 
survey to the bubble imprint. 

Figures [2] [3] and [4] show contour plots for models A, B, and 
C, denoting regions of parameter space where our surveys are able 



0.15 







0.05 



bub 



60 
(Mpc) 



100 



Figure 3. Model B: Contour map of detection. As for Figure [2] For com- 
parison with other figures, note that G(z = 6)/G(z = 0) = 0.18, so that 
[G(z = 0)/G(z = 6)]e b lies in the range [0,1]. 




bub 



60 
(Mpc) 



100 



Figure 4. Model C: Contour map of detection. As for Figure [2] Note the 
greatly decreased ability of SDSS alone to detect bubbles when compared 
with Figure|2] 



to make a detection. We shade the region of the (rbub, e&) plane 
where a detection can be made by SDSS alone (white), SDSS and 
S2 combined (grey), and where no detection can be made (black). 
For models A and C, we consider the region of parameter space 
with r bu b € [5 Mpc, 100 Mpc] and e b G [0.1,1]. For model B, we 
choose a slightly different normalisation so that [G(z — 0)/G(z — 
6)]e(, £ [0.1, 1]. This makes the range of amplitude of Pbub at 
z = identical in the region covered in Figures [2] and [3] Note that 
G(z = 6)/G(z = 0) = 0.18 in our fiducial cosmology. 

First, compare Figures|2]and|4] In model A, we see that SDSS 
alone is able to detect bubbles over a wide range of tb, provided 
that the bubbles are large (rbub > 40 Mpc). In contrast, when we 
allow the bubble amplitude to decrease with time, as in Model C, 
we see that SDSS alone is almost unable to constrain either bub- 
ble parameter. In both cases, addition of the S2 survey greatly im- 
proves the situation, allowing a wider range of parameter space to 
be probed. However, even with S2, the theoretically preferred re- 
gion with rbub < 10 Mpc and et < 0.2, towards the bottom left 
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hand corner, remains unconstrained. The prospects for detection are 
clearly enhanced by including galaxy surveys at higher redshift. 

Figure [3] for Model B, shows that growth improves the 
prospects for probing smaller values of tb, but makes little differ- 
ence to our ability to constrain the bubble size. As in model A, 
SDSS alone can only probe bubbles with rbub > 40 Mpc, and S2 
is required to probe smaller scales. Note that when we normalize 
to the present day, the inclusion of growth in model B reduces the 
amplitude of the bubble imprint seen by the S2 survey by a factor 
of [G(z = 3)/G(z = 0)] 2 « 0.1 over that in model A. This is re- 
sponsible for the increased region that is undetectable to SDSS+S2 
in Figure [3] The amplitude of Pbub at z = 0.3 is very similar 
in these two models, resulting in the nearly identical contours for 
SDSS only (when rescaled to account for the different normalisa- 
tion) in Figure s|2] and [3] The striking differences between the three 
models indicates the importance of the time evolution of the bubble 
imprint. 

It is worth pausing for a moment to consider where our lever- 
age on the bubble power spectrum originates. When we combine 
the two surveys, the bulk of the improvement is coming from the 
S2 survey alone. This is unsurprising, as the growth of the density 
fluctuations means the bubble imprint is a more significant contri- 
bution to the galaxy power spectrum at early times. Further, if we 
consider FigureQ] we see that for very small bubble sizes, the bub- 
ble spectrum begins to resemble white noise over the region probed 
by the galaxy surveys. This would further complicate detecting the 
bubble imprint as it could then be confused with residual Poisson 
shot-noise in the galaxy counts. This problem is greatest at low z 
where the non-linear scale is larger. Both of these motivate per- 
forming this test in galaxy surveys at increasing redshift, ideally at 
the redshift of reionization, where an Ha survey may be possible. 

We conclude that galaxy surveys should be sensitive to the 
imprint in the galaxy power spectrum left over after reionization. 
However, detecting this imprint will be difficult unless the charac- 
teristic size of HII regions is large (rbub > 10 Mpc) and the effects 
of feedback significant (e;, > 0.1). This should be sufficient to con- 
strain the more extreme models for reionization, but is unlikely to 
impact more reasonable scenarios. There is significant uncertainty 
in this prediction stemming from the difficulty in predicting the 
evolution of the imprint to more recent times. 

Currently, the best hope for measuring the size of HII 
regions during the early stages of reionization lies with up- 
coming 21cm observations (e.g., LOFAeQ MW^E or PASTO). 
Direct imaging of the HII regions is unlikely with the first 
generation of detectors, but the prospects for statistical detec- 
tion at z < 10 are good (Zaldarriaga, Furlanetto & Hernquist 
l2004l ; lBowman. Morales &~ Hewitt 2006; McOui nn et alj2005h . At 
higher redshifts, z > 10, corresponding to lower frequencies, sky 
noise increases dramatically making observations more difficult. 
An imprint directly upon the galaxy power spectrum avoids these 
technical issues, making possible a complementary measurement. 
In the event of very early reionization, detection of the imprint dis- 
cussed in this paper might provide weak constraints on reionization 
before 21cm experiments reach the desired sensitivity. 



1 See http://www.lofar.org/ 

2 See http://web.haystack.mit.edu/aiTays/MWA/ 



See Pen. Wu & Peterson 1 2005) 



6 IMPLICATIONS FOR DARK-ENERGY CONSTRAINTS 

Current constraints on dark-energy parameters arise from the com- 
bination of high-precision CMB data with information from galaxy 
surveys. The combination of high-jz (z > 3) information, long 
before dark energy becomes dynamically important, with low- 
z (z < 3) information, deep within the dark-energy-dominated 
regime, serves to break many of the degeneracies that either data set 
possesses when used alone. Adding in more galaxy surveys at dif- 
ferent redshifts further constrains the evolution of the dark energy, 
allowing constraints on both fix and its equation-of-state parame- 
ter w(z). In the previous section, we considered the bubble imprint 
as a useful signal; in this section we consider it as a potential source 
of noise for galaxy surveys. If the bubble power spectrum is able to 
mimic the effects of dark energy, then it will degrade our ability to 
constrain dark-energy parameters. Throughout this section, we will 
consider a dark-energy model with wo = — 1 and w\ = 0. Our 
numerical results depend upon this choice of model, but the overall 
picture remains the same when wo and Wi take other values. 

Galaxy surveys provide direct constraints on the dark en- 
ergy through both the baryon oscillations and from the Alcock- 
Paczynski effect. They also provide indirect constraints in combi- 
nation with CMB data, as they probe fl m independently of O m /i 2 , 
the parameter directly probed by the CMB. This allows the CMB 
indication of flatness Qk ~ to constrain fix- The bubble imprint 
must interfere with one of these measurements to be a source of 
confusion. 

Measurement of the baryon oscillations allow a determination 
of the angular-diameter distance. Their distinctive oscillatory struc- 
ture is very different from the smooth structure that we expect from 
any plausible bubble imprint and so we do not expect there to be 
any confusion between the two. The inferred peak position, ampli- 
tude, and overall shape of the galaxy power spectrum, on the other 
hand, could be affected by the smooth form of the bubble imprint, 
making these the most likely points of confusion. Thus, we expect 
parameters such as fl m and n s to be sensitive to the bubble imprint. 
This simple picture is modified by inclusion of CMB data, which 
places tight constraints on many of these parameters making the 
effect of the bubble imprint more subtle. 

The correlation between the different parameters is indicated 
in Figure|5] for a model with rbub = 10 Mpc and et = 0.1. Note 
that there is a weak correlation between rbub and the dark-energy 
parameters. A slightly larger covariance is seen between rbub and 
n s . At larger values of (rbub, Cb), the picture remains unchanged 
except for a breaking of the degeneracy between rbub and eb as the 
cutoff in -Pbub(fc) on small scales becomes more pronounced. 

Before detailing the effect the bubble imprint has on statistical 
errors, let us consider the possibility of systematic biasing of our 
best-fitting values, if an existing bubble imprint was ignored in the 
analysis of data. This will be relevant only in the case that the bub- 
ble imprint is not easily detectable, as an obvious imprint would 
certainly be included in the data analysis. For the case where the 
bubbles are not detected - i.e., e;, < 1a^ h and rbub < 2oy bub - we 
have estimated this systematic offset between the inferred and true 
parameters, using the Fisher matrix to approximate the full likeli- 
hood surface. We find that the offset is significantly smaller than the 
parameter uncertainty, typically being of order ~ 0.1%. From this, 
we conclude that failing to include the imprint should not system- 
atically affect parameter estimates in the near future. When galaxy 
surveys begin to probe cosmological parameters below the percent 
level this effect will need to be included. We now turn to the effect 
of the imprint on parameter constraints. 
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Figure 6. Model A: Contour map of errors in wo in the bubble parameter 
plane. We plot contours spanning the range ff WQ = 0.22-0.58 in intervals 
of 0.04. The fiducial model takes too = — 1 an d w i = 0- When no bubbles 
are present, we find cr WQ = 0.39. 



U m h 2 Ux wo toi n s A 2 S r bub e b 



Figure 5. Illustration of the reduced covariance matrix. 
(Q m /i 2 ,&x,WQ,wi,n s ,A^,rb\ibi e b)- The model uses rb u b = lOMpc 
and ef, = 0.1. Black indicates strong correlation and white indicates little 
correlation. 



Figures [6] [7j and [8] indicate error contours for wo over the 
(Ttmb, £fc) plane. The same shading scheme is used in all three fig- 
ures to allow easy comparison. First consider Figure[6] We see that 
the uncertainty on wo is maximal for bubble parameters rbub ~ 
80 Mpc and ej, ~ 0.5. The form of Pbub{k) is plotted in FigureQ] 
where we see that it cuts off close to the maximum of the density 
power spectrum. This is consistent with our above statements. We 
find a maximum uncertainty of a WQ — 0.48 in contrast with the 
uncertainty a wo = 0.39 in the absence of bubbles. This indicates 
that bubbles can be an important source of noise in attempts to con- 
strain dark energy. However, the large values of (rbub, £b) required 
for this effect seem theoretically unlikely and from the discussion 
in SJ5] would allow direct detection of the bubbles. For more rea- 
sonable choices of bubble parameters (rbub < 10Mpc,£b < 0.1), 
the uncertainty on u>o reduces to a WQ = 0.39. Thus, the effect of 
the bubble imprint is likely to be somewhat important in future at- 
tempts to constrain dark energy. 

Now consider Figures [7] and [8] The increased uncertainty in 
wo is more pronounced in Model B, where the uncertainty rises as 
high as a wo = 0.62. Even here, for small values of (rbub, £(>), the 
bubble imprint becomes unimportant and we recover a WQ — 0.39, 
the no-bubble uncertainty. Again this maximal uncertainty occurs 
close to rbub ~ 80 Mpc and e;, ~ 0.5. The increased value of a wo 
is a consequence of the bubble imprint growing at the same rate 
as the density fluctuations. Consequently, the overall shape of the 
total galaxy power spectrum remains constant in time. Thus, com- 
bining information at two redshifts provides much less leverage on 
separating out the two components, leading to larger parameter un- 
certainties. In model C, the damping of the imprint means that it 
has much less effect on the dark-energy parameters. 

Finally, we note that Figures [6] |7J and [8] display a region of 
decreased uncertainty in wq in the top right hand corner, when the 
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Figure 7. Model B: Contour map of errors in wo in the bubble parameter 
plane. As for Figure [6] 



bubbles are large and feedback strong. This is an interesting ex- 
ample of the AP effect. In this region, the bubble power spectrum 
dominates over the density contribution and the overall shape of the 
galaxy power spectrum displays a well defined, sharp cutoff. Dis- 
tortion of this scale by the AP effect places good constraints on the 
dark energy. Galaxy surveys already show that the galaxy power 
spectrum closely traces the underlying density field, so this region 
is ruled out. 

Having considered dark-energy parameters, it would seem 
natural to also consider inflationary parameters; e.g., the tilt n s 
and amplitude A%. For the surveys that we have analysed, inclu- 
sion of the bubble power spectrum makes little difference to the 
uncertainty on these parameters. Essentially, all of the informa- 
tion needed for constraining these quantities is contained within 
the CMB. In the absence of information on the optical depth r, 
or if there are significant tensor modes, galaxy-survey information 
becomes important in breaking degeneracies. This is not true in 
the cases that we consider, where CMB -polarization information 
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Figure 8. Model C: Contour map of errors in too in the bubble parameter 
plane. As for Figure [6] 



is well measured. If we were to try and use galaxy-survey data by 
itself, we would notice increased uncertainty in the tilt n s . 



7 CONCLUSIONS 

In this paper, we have discussed the possibility that patchy reion- 
ization may leave an imprint in the distribution of galaxies through 
its effect on the collapse and cooling of baryons. We considered a 
simple ansatz linking galaxy number density to the ionization frac- 
tion and used a halo-model approach to calculate the imprint of 
inhomogenous ionization on the galaxy power spectrum. We then 
applied a Fisher-matrix approach to place constraints on the effect 
of this imprint. 

Our calculation shows that detecting the bubble imprint 
through large galaxy surveys is potentially feasible, but highly de- 
pendent upon the details of reionization. We have shown that, for a 
detection to be possible with upcoming experiments, bubbles must 
be large (r"bub > 10 Mpc) and the feedback moderately strong 
(eb > 0.1). This suggests that the most reasonable region of pa- 
rameter space (rb u b < 10 Mpc, tb < 0.1) will not be detected 
with currently proposed galaxy surveys at z < 3. Potentially, a 
z ~ 6 galaxy survey might give the additional leverage needed for 
a concrete detection. 

Beyond the possibility of detection, we have considered the ef- 
fect of the bubble imprint on constraining dark-energy parameters. 
We find that the distinctive nature of the baryon oscillations helps 
minimize any degeneracy arising. Only if the characteristic bub- 
ble size is ~ 80 Mpc does the bubble imprint seriously impact our 
uncertainty in wo. In this case, the bubble power spectrum closely 
mimicks the cutoff of the density power spectrum. This is a region 
of parameter space where the bubbles should be easily detected. 
For more sensible values of (rb u b,£t), there is little or no impact 
on dark-energy constraints. When the bubbles are not detectable, 
we find that ignoring them in the analysis of galaxy data does not 
introduce any significant biasing of the best-fitting parameters. 

Our approach has emphasised the use of a simple toy model 
to probe the effect of reionization on the distribution of galaxies. 
If future galaxy surveys are able to make detections of this signal, 
it will be important to incorporate more detailed physics to better 



constrain the shape and amplitude of the bubble imprint. With our 
present understanding of reionization, this seems premature. 

Future galaxy surveys will greatly add to our knowledge of 
the distribution of galaxies and the nature of the dark energy. If we 
are to extract maximum information from these surveys, we must 
tighten our understanding of the biasing of galaxy formation and 
the possible effect of reionization on early generations of galaxies. 

This work was supported at Caltech in part by DoE DE-FG03- 
92-ER40701 and NASA NNG05GF69G. 
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